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ABSTRACT 

This  study  considered  the  feasibility  of  obtaining 
a  self-adapting  pitch  rate  control  system  of  a  hypothetical 
airframe,  utilizing  linear  compensation  of  the  invariant, 
i.e»,  normal  operating  state,  parameters  of  the  aircraft. 

The  relationship  between  the  flight  dynamics  and 
aircraft  parameters  was  established,  and  each  was  reduced 
to  an  invariant  and  variant  component.  A  method  of  repre- 
senting the  system  with,  and  without  parameter  variation 
was  then  developed.  Linear  compensation  of  the  normal  op- 
erating state  was  developed  to  provide  nearly  flat,  minimum 
phase  shift,  frequency  response;  as  the  primary  performance 
criterion. 
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The  compensated  and  uncompensated  aircraft  was 
analyzed  in  widely  varying  operating  conditions.  Theoret- 
ical frequency  response  and  analog  computer  transient 
response  simulation,  yielded  the  following  results: 

The  desired  frequency  response  in  all  flight  con- 
ditions can  almost  be  obtained,  since  the  compensated 
system  introduced  varying  attenuation  with  a  maximum 
of  35%.      Six  of  the  seven  conditions  produced  acceptable 
transient  response. 

To  eliminate  the  attenuation,  a  variable  gain 
amplifier,  sensitive  to  the  conditions  is  proposed.  Us- 
ing linear  compensation  the  complete  system,  including 
the  model,  amplifier,  and  pitch  rate  system,  acceptably 
reproduced  the  desirable  model  output,  through  most  of 
the  various  flight  conditions.   This  satisfies  the 
basic  requirements  of  a  self-adaptive  system. 

Thesis  Supervisor:   Y..T*  Li. 

Title:   Associate  Professor  of  Aeronautics  and 
Astronautics. 
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OBJECT 

The  object  of  this  study  is  to  propose  and 
investigate  a  self-adaptive  control  system  by 
means  of  compensation  of  normal  operating  state 
aircraft  parameters,  using  linear  techniques. 


CHAPTER  1 
INTRODUCTION 


Modern  Control  System  Requirements 

With  the  advent  of  high  performance  aircraft* 
satellites,  rocket  aircraft,  guided  missiles,  ICBM's, 
and  eventually,  space  ships,  there  has  been  an  acceler- 
ated requirement  for  replacement  of  the  human  pilot  with 
automatic  flight  control  systems.   This  replacement  is 
not  necessarily  complete,  and  in  many  cases  is  used  only 
to  augment  the  human  characteristics  in  order  to  add  damp' 
ing, increase  accuracies,  and  in  extreme  cases,  to  provide 
control  capability.   This  class  of  automatic  control 
system  has  been  further  instrumented,  and  is  more  popu- 
larly known  as  a  self-adapting  control,  or  servo,  system. 

The  human  servo  system  is  an  excellent  example  of 
self-adaptation.  However,  in  the  present  day  supersonic 
vehicles,  reaction  time  is  too  slow,  and  muscular  coordin- 
ation is  not  sufficient  to  cope  with  the  various  extremes 
of  velocities,  temperatures,  forces  and  other  modifying 
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influences;   since,   modern  aircraft  have  flight  capa- 
bility over  extreme  ranges  of  altitude,  air  speed  and 
Mach  number.   This  capability  is  mostly  attainable  due 
to  propulsive  thrusts  which  approach  maximum  aircraft 
gross  weight.   Automatic  control  systems  and  stability, 
augmentation  devices  have  thus  been  installed  in  aircraft 
to  provide  satisfactory  stability  and  control  throughout 
the  various  flight  regimes.   The  stability  and  control 
characteristics  of  all  aircraft  vary  as  functions  of 
equivalent  air  speed,  Mach  number,  altitude,  aircraft 

gross  weight,  center  of  gravity  position  and  external 

(1) 
configuration. 

System  Representation  and  Classification. 

Self-adaptation  may  often  be  represented  by  an 
ordinary  negative  feedback  control  loop  which  in  itself, 
is  the  first  of  the  basic  servo  concepts.   However,  the 
quality  of  this  adaptation  may  then  be  used  as  a  measure- 
ment of  performance  capability. 

Some  basic  classifications  of  the  adaptive  concept: 
The  first,  and  most  easily  understood  is  that  of  negative 
feedback.   If  there  is  a  high  gain  in  the  forward  loop, 
assuming  non  unity  feedback  then  the  closed  loop  transfer 
function  then  becomes  approximately  that  of  the  negative 
feedback  performance  equation.   The  forward  loop  has  no 


effect  essentially,  on  the  entire  closed  loop  if  the 
combined  gains  are  much  greater  than  one.   Thus,  if 
the  forward  loop  is  the  controlled  member,  then  the 
system  operation  will  be  completely  independent  of 
the  controlled  member,  e.g.,   an  aircraft,  which  is 
an  ideal  situation  in  self-adaptation.   However, 
actual  implementation  with  mechanical  or  electrical 
devices  are  not  as  easily  found  as  the  transfer 
equation  would  have  us  believe. 

There  are  many  ways  of  re-writing  the  block 

diagram  representing  the  above  concept  in  order  to 

overcome  the  basic  difficulties  with  normal  feedback. 

(2,3) 
The  method  that  is  introduced  by  a  number  of  authors  . 

is  that  of  the  conditional  feedback  loop  based  on  a 
performance  model,  that  is  designed  such  that  the  over- 
all system  approximates  the  characteristics  of  the 
idealized  model,  over  wide  ranges  of  operation. 

(2) 

Sperry  Gyroscope  Company  have  further  described 

systems  such  as  conditional  feedback  with  a  pre-filter, 
which  will  be  further  discussed  in  the  later  chapters 
of  this  reporto 

Compensation  and  Effects. 

Automatic  control  systems  are  generally  multiple 
closed  loop  with  feedback  paths  utilizing  variously: 


displacement,  velocity,  and  acceleration  of  the  aircraft 
as  control  quantities.   Therefore,  the  static  and  dynamic 
response  characteristics  have  a  great  effect  and  cannot 
be  ignored  by  using  only  high  'gain  feedback,  since  satura- 
tion and  power  obtainable  from  the  physical  system  often 
falls  far  short  of  the  theoretical  desires.   It  becomes 
obvious,  then,  that  some  compensation  must  be  made  for 
the  aforementioned  effects. 

This  compensation  may  take  many  forms  as  predi- 

(5) 
cated  by  Whitaker,  Yamron,  and  Kezer; .  namely,  changing 

system  sensitivities  as  calibrated,  programmed  functions 
of  indicated  airspeed,  Mach  number  and  other  flight  quanti- 
ties.  The  adjustments  to  the  system  sensitivities  are 
then  made  such  that  they  are  open  loop,  from  the  stand- 
point of  the  effect  of  the  adjustment  on  system  perform- 
ance. Past  action  adapting  control  systems  are  now  be- 
coming of  greater  importance,  as  stated  before,  and  will 
become  even  more  so  with  the  advent  of  space  flight. 
There  will  be  colossal  changes  in  environmental  parameters 
as  well  as  the  basic  system  parameters  themselves.  Large 
changes  in  mass  with  attendant  changes  of  inertia,  aero- 
dynamic heating  due  to  high  velocities  in  the  earth's 
atmosphere,  and  changes  of  the  density  of  the  medium  in 
which  the  controlled  vehicle  is  operating,  are  a  few 
examples  of  the  changes  that  must  be  anticipated  and 
corrected  for  by  the  vehicles  control  system.   Recently 


it  was  announced  that  the  problem  of  re-entry  of  a 

manned  satellite  would  be  undertaken.   Let  us  assume 

that  the  satellite  would  re-enter  at  an  altitude  of 

100  kilometers,  and  utilize  a  bounce-glide  technique 

of  re-entry,  i.e.,  decelerate  gradually  in  the  upper 

layers  of  the  sensible  atmosphere  until  a  tolerance  of 

temperature  is  obtained,  then  bouncing  back  to  the  next 

higher  level  and  then  descending  again,  in  gradual 

sweeps,  for  re-entry,.   This  requires  that  a  theoretical 

aircraft,  weighing  20  tons  with  a  wing  loading  of  30 

kilograms  per  square  meter,  undergoes  skin  temperature 

variations  of  from  800  degrees  C,  to  400  degrees  C, 

to  1400  degrees  G,  to  0  degrees  C,  while  undergoing  a 

velocity  change  from  8.5  kilometers  per  second  to  zero 

(4) 
during  a  2-hour  flight. 

It  is  seen,  then,  that  it  is  necessary  either  to 

change  certain  system  parameters;   or  to  modify  command 

signals,  as  the  environment  changes,  in  order  to  maintain 

satisfactory  dynamics.   It  is  further  necessary,  since 

flight  testing  in  outer  regions  of  space  vehicles  is  not 

physically  realizable,  to  be  able  to  design  a  theoretical 

system  that  can  be  implemented  at  a  later  date  that  is 

intelligent  enough  to  know  when  dynamic  performance  has 

changed  and  to  utilize  this  knowledge  to  a  dapt  itself, 

so  that  its  final  performance  will  be  within  a  satisfact- 

(3) 
ory  range  at  all  times. 


The  present  trend  to  supplement,  or  even  eliminate, 
a  human  pilot,  has  led  to  severe  requirements  for  control 
performance,  since  the  capability  of  many  applications  re- 
quires control  through  violent  maneuvers  and  yet  must  still 

(10) 
be  sensitive  to  small  input  command  signals.    In  many 

cases  a  control  system  may  be  extremely  complex;   however, 

when  dealing  with  a  specific  vehicle  such  as  an  aircraft, 

we  are  limited  by  size,  weight  and  structural  consider- 

ationso   The  result,  then,  is  often  a  servo  with  a  limited 

output  rate.   If  we  are  able  to  make  this  assumption,  the 

controller  may  be  considered  as  a  simple  linear  system 

(14) 
with  a  limit  on  the  rate  of  control.    But,   it  is  this 

very  limit  on  the  control  surface  rate  that  is  a  major 
difficulty  in  a  system  requirement  for  stability  for  both 
large  and  small  commands.   Compensation  and  adjustment 
of  the  system  parameters  to  give  the  required  stable  res- 
ponse generally  results  in  undesirable  slow  response  to 
small  signals.   There  have  been  many  theories  on  compen- 
sation for  such  undesirable  effects,  more  specifically: 
to  design  system  parameters  as  functions  of  error  and 
error  rates.   This  has  led  to  development  of  non-linear 

elements  to  improve  the  performance  of  both  linear  and 

(5) 
non-linear  devices.-  . 

There  have  been  methods  proposed  to  also  accomp- 
lish this,  in  which  a  control  system  tends  to  adjust 


itself  to  a  pre-determined  operating  condition..  The 
system  imposes  known  continuous  variations  on  system 
inputs  to  determine  system  parameters,  to  maintain  out- 
put at  a  maximum. 

Series  compensation  and  unity  feedback  with  a 
controlled  component  in  which  the  series  compensation 

is  varied  by  an  external  computer  is  yet  another  propo- 
(6) 

sal. 

Summary 

Thus,  adaptive  mechanisms  may  be  cataloged  into 
two  general  types,  the  first  of  which  is  the  simple  high 
gain  system  using  linear  negative  feedback,,  The  second, 
which  is  of  primary  importance  in  the  recent  history  of 
self-adaptation,  is  where  the  alteration  of  the  control- 
ler as  a  function  of  system  response,  is  involved  in  a 
non-linear  system.  This  type  leads  ultimately  to  an  op- 
portunity to  design  and  manufacture  a  controller  that 
is  optimum  insofar  as  the  continuous  determination  of 
the  static  and  dynamic  characteristics  of  the  controlled 
members  is  concerned. 

To  correct  the  difficulties  in  stability  already 

mentioned,  it  has  been  s uggested  that  uniform  stability 

may  be  obtained  by  compensating  for  changes  in  aircraft 

(12) 
stability  directly,,     This  is  accomplished  by  adjusting 
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control  parameters  as  functions  of  response  in  order 

(8) 
to  maintain  a  uniform  response.   The  basic  instabilities 

generally  result  from  high  feedback  gain,  due  to  the 
dynamic  effects  of  the  servo  actuators.  Todays  manufact- 
uring techniques,  however,  have  been  able  to  synthesize 
controllers  with  better  rates  of  response,  and  thus 
higher  gains  are  obtainable,  leading  to  more  practical 
value  in  the  theory  of  negative  feedback. 

By  definition,  optimum  control  of  the  aircraft 
results  if  the  control  system  provides  close  control  of 
the  transient  and  steady  state  responses  of  the  aircraft 
to  command  and  disturbance  inputs,  so  as  to  use  the 
capabilities  of  the  aircraft  most  effectively  for  ful- 
filling the  specifications  of   its  flight  mission.  The 
characteristics  involved  being  response  time,  damping, 

dynamic  and  static  errors,  and  control  of  interference 

(2) 
effects.   This  very  definition  allows  great  flexibility 

in  utilizing  the  aircrafts  capability,  without  forcing 
an  impractical  response  from  the  physical  vehicle.  There 
are  as  many  design  criteria  as  there  are  designers,  how- 
ever, and  in  selecting  such  a  criteria  it  must  be  remem- 
bered that  the  desired  performance  itself  may  vary  from 
one  phase  of  a  flight  to  another. 


CHAPTER  2 
SCOPE 

The  purpose  of  this  study  is  to  investigate  the 
feasibility  of  obtaining  a  self-adaptive  system  by  em- 
ploying a  linear  technique,  herein  defined  as  invariant 
parameter  compensation.   The  basic  philosophy  underlying 
this  technique  is  the  establishment  of  a  state  in  which 
the  airframe  parameters,  and  hence  the  airframe  dynamics 
are  defined  to  be  invariant.   .Further,  having  established 
the  airframe  in  this  state,  it  is  then  assumed  that  suit- 
able compensation  can  be  devised  to  produce  an  essentially 
flat  response  within  specified  tolerances,  of  the  air- 
frame controlled  system  over  a  wide  range  of  frequencies. 

The  design  tolerance  criteria  is  as  follows? 

a.  Amplitude  ratio  l-jfcljfc  variation 

b.  Phase  shift  of  less  than ±10° 

all  within  an  excitation  frequency  range  of  0-50 
radians/second. 

The  value  of  this  technique  lies  in  the  capability 
of  establishing  an  effective  natural  frequency  very  much 
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greater  than  any  frequency  to  which  the  system  will 
normally  be  subjected,  and  an  effective  damping  ratio 
which  is  dynamically  acceptable.   Through  this  technique, 
it  is  believed  that  the  normal  range  of  parameter  varia- 
tions to  which  the  aircraft  will  be  subjected  will  not 
change  the  performance  function  to  any  great  extent, 
over  the  range  of  interest  of  input  frequencies.   It 
is  further  desirable  that  the  aircraft  perform  in  a 
manner  as  specified  by  the  performance  of  a  reference 
model  system  which  precedes  the  controlled  loop,  and 
whose  output  will  be  the  input  to  the  controlled  loop0 
The  performance  function  of  the  model  will  be  developed 
in  later  chapters. 

If  the  aircraft,  as  compensated,  does  reproduce 
as  its  output,  the  output  of  the  reference  model  system 
for  all  normal  variations  of  parameters,  and  for  all 
frequencies  within  a  range  of  interest;  then  the  aircraft 
control  system  will  truly  be  self-adaptive. 

The  scope  of  this  paper  will  be  limited  to  the 

pitch  rate  performance  of  a  hypothetical  aircraft,  whose 

performance  function  and  dynamics  are: 

Performance  function  Diagram: 

AIRFRAME 


a 
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Whe re : 

PP  _  Reference  model  system  performance  function, 

PFS^A=  Servo  and  Actuator  performance  function,, 

Za =  First  order  zero;  an   airframe  dynamic. 

Wn=  Airframe  natural  frequency;  an   airframe 

dynamic. 

$   =  Airframe  damping  ratio,  an   airframe 

dynamic o 

0 

9  in  —  Commanded  pitch  rate. 
0  m- Model  commanded  pitch  rate, 
o  E=Elevator  deflection  angle. 
0O  ^Airframe  pitch  rate  response. 

In  order  to  employ  this  method  of  obtaining  self- 
adaptation,  it  will  be  necessary  to  establish  the  relation- 
ship between  the  above  airframe  dynamics  and  their  correspond- 
ing parameterso  Further,  since  compensation  is  to  be  performed 
on  only  the  invariant  parameters,  it  is  necessary  that  a 
method  be  devised  wherein  the  invariant  parameters  may  be 
handled  separately  from  the  variant  components. 

This  development  will  lead  to  representation  of 
parameter  variations  by  synthesis  of  a  fictitious  system 
containing  both  the  invariant  and  variant  components.  The 
variant  components  represented  in  this  fictitious  system 
are  the  inner  loops,  which  are  active  only  when  the  air- 
craft is  not  in  its  invariant  state. 
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Although  it  is  anticipated  that  this  technique 
can  be  utilized  by  any  aircraft  configuration,  a  hypoth- 
etical airframe  was  employed,  since  by  this  method  the 
contributions  o£   the  various  parameters  to  the  aircraft 
dynamics  and  the  effect  of  variations  of  these  parameters 
on  the  dynamics  could  be  more  adequately  demonstrated. 

The  invariant  system  parameters  are  those  which 
exist  when  the  subject  aircraft  is  in  the  state  of  alti- 
tude, airspeed, . and  loading  specified  as  the  invariant 
state  of  operation.  The  specification  of  this  state  appears 
to  be.  somewhat  arbitrary,  although  it  would  seem  logical 
to  designate  the  "normal"  operating  state  of  the  aircraft 
as  its  invariant  state,  l"br  aifrarnes  not  possessing  such 
a  state,  selection  of  the  invariant  parameters  becomes 
a  matter  of  additional  study  and  specifications  of  criteria 
to  establish  the  "best"  set  of  invariant  parameters  for 
such  types  of  vehicles  as  ballistic  missiles,  may  well 
prove  to  be  of  sufficient  magnitude  to  justify  further 
work  in  this  field. 
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CHAPTER  3 


THE  HYPOTHETICAL  AIRFRAME 


Definition  of  Parameters 

As  indicated  in  Chapter  2,  the  effectiveness  of 
linear  compensation  is  to  be  tested  on  a  hypothetical  air- 
frame.  In  this  chapter  the  airframe  pitch  rate  equation, 
dynamics  and  parameters  are  established.  In  addition,  the 
invariant  state  will  also  be  defined. 

For  the  purpose  of  this  development,   the  following 
table  of  parameters  is  applicable: 

Dimensional 

f>  —  air  density 

f>  —  air  density  at  sea  level 

V  —  airframe  airspeed 

S  =  a  characteristic  surface  area  of  airframe 

c  —  a  characteristic  chord  length  of  airframe 

Orr\  —  moment  coefficient 

Ct  =  coefficient  of  lift 

C«  =  derivative  coefficient  =dCi  ;  i  :=m,L  ;  j-^,*,^* 
J  dj 
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£^  =  control  surface  deflection 

c<  =  airframe  angle  of  attack 

Q     =   airframe  pitch  angle 

Q   =  airframe  pitch  rate;  i.e.  d© 

dt 

2 
Iyy=  longitudinal  axis  moment  of  inertia  =  km 

k   =  longitudinal  axis  radius  of  gyration 

m   =  mass  of  airframe 

Non-dimensional 


r  =jr    >    *  =^    '     «>  =  * 


Aircraft  Performance  function 

The  short  period  performance  function  relating 
elevator  angle  (6^)  input  and  aircraft  pitch  rate  (60) 
output  may  be  represented  in  terms  of  dynamics  as  follows: 

4         -  kI>*aQ 


4  .  \f+ZT*M>  *  <] 

(Note:   The  negative  relationship  is  introduced 
by  sign  conventions  for  Se  9  0o     •         ^  should  be  apparent 
that  the  same  sign  convention  will  apply  for  Q~tQi    ,    thus 
making  a  positive  relationship  for  Q?   r    Bo  * ) 
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It  is  assumed  that  each  of  the  dynamics  my  be  repre- 
sented as  the  sum  of  a  variant  and  an  invariant  component, 
as  follows: 

Za  =  Zao  +  A  Za  ;   K  =  Ko  4-  4  K  ;  /Wn  =  (JWn)0+  4  (/Wn) 
wn=  {Wn}o  +  A(wn> 

The  subscript  zero  signifies  the  parameter  is  in- 
variant.  The  selection  of  Wj?  &  3*Wn  as  basic  parameters 
was  done  to  eliminate  cross  product  terras  in  later  develop- 
ment. 

Dynamics  in  terms  of  airframe  /flight  path  parameters 

w*     =-AVV*Acry     where.     Al  =  ^f-,    ;A2  =  £*Sl£* 

2  .f  Wn  -  -A      \J~  where  ;      A,  =  -  [C">Q  *  Cm*  / 1*^* 

^  3  L  ic  K*  A 

Za    =     A  ^  where;      A4  =  glft^  C^e  -  G*£<1 

=   ^—\?'A^  where;      A7=_.  =   ^^ 

^  2c^;'  8     *c^€ 

The  following  parameters  are  assumed  to  be  constant: 

2 
S>  c,  K„,  Gyy,^   ,  C^^  ,  C/.^  ,  C^,^  ,  £>Ls£     y    &  Cm^p  <•   Actually, 

the  derivative  coefficients  vary  somewhat  with  Mach  number. 
However,   the  variations  of  these  quantities  are  not  signifi- 
cant when  compared  to  variations  of  other  parameters.  Conse- 
quently,  the  results  obtained  will  not  be  too  great  in 
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error  as  a  result  of  this  simplifying  assumption. 

Just  as  the  airframe  dynamics  are  considered  to 
be  composed  of  a  variant  ana  an  invariant  component,  so 
too,  may  the  airframe/f light  path  parameters  be  considered, 
The  following  then  define  the  dynamic  components  in  terms 
of  the  corresponding  parameter  components: 


(a)  ul*   =  (u>£)    +   4  a/* 


o 
2 


(b)  */">*  =  2f/^jQ  +  Z^(fuj^) 

(c)  f*.  =  F«  *  42?«, 

ac    =  <,  *   4  a; 
^  /T     =  , — * ^  (AFC  V0  +A\Z)\<Z  +  *<r)(."o  +4S*)'  Vc  VA\ 
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Specification  of  the  Invariant  State 

The  hypothetical  airframe  is  assumed  to  be  at  its 
invariant  state,  i.e.  normal  operating  state,  at  an  alti- 
tude of  35,000  ft.,  with  an  airspeed  of  1000  ft. /sec.  and 
a  mass  of  497  slugs.   This  condition  produces  the  following 
values,  defined  to  be  the  invariant  values,  of  the  parameters: 

fi   _  .835  X  103  slug/ft3   <r^   =    .35 

Vo  =  lOOOft/sec.  u^      =         79.1  (see  below) 

Mo  =  497  slugs 

In  addition,  the  parameters  previously  defined  as 
constants  have  been  assigned  the  following  values: 

Parameters  Coefficients 

k|=  81  ft2  Ax  =  -0.00275 

k|=  1.046  A2  =  0.0485 

£  =  300  ft2  A3  =  -0.475 

c  =  8.8  ft  A4  =  -0.486 

C,^=  -4.5  A5  =  -2.0 


<->»«•<  • 


-.45  A«  =  -0.6 


CU=   -3.0  A?  =   0.00306 

0*9=   -7.0  Ao  =  -0.30 


Ct?e=   -0.2 
0^=   -0.5 
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Utilizing  the  foregoing  parameters,   the  invariant 
dynamics  are  determined  to  be: 

(W|)0=  13.1 

2(fWn)0=  2.1 

(Za)Q  =  1.07 

(K)Q   =  15.55 

Once  the  dynamics  have  been  defined,  the  problem 
is  reduced  to  one  of  compensation;  to  obtain  acceptable 
aircraft  performance. 
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CHAPTER  4 

COMPENSATION  AND  EFFECT 
OF  PARAMETER  VARIATIONS 


Method 

As  previously  indicated,  this  study  is  concerned  with 
compensating  the  air  frame  for  one  state  of  operation  and 
then  testing  the  effectiveness  of  this  system  when  the  vehicle 
is  in  states  other  than  that  in  which  it  is  compensated. 
It  was  found,  that  handling  of  variational  effects  could  be 
greatly  facilitated  by  employing  the  analog  representation 
of  variant  and  invariant  parameters  shown  below. 

(a)  Functional  Diagram  and  Exact  Equation  of  -Sa.  : 
Diagram  of  elevator  to  aircraft 


Define  KG-A/C  such  that: 

(p2+  2jWn  p  +  W2)   qQ  =  K  (p  -+-  Za)<f£ 

(1)   [p2+  2([VJ.+  *M  )P+(W§)0+A(W|)]  q0  ^ 

(Ko+AK)  (p-4- Zao+aZa)  o~e 

Where  the  subscript  o  refers  to  the  invariant 

parameters 
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(b)  Analog  Representation  of  Equation   (1) 


j£> 


fSx 


<ksi    H$ 


"^ 


^ 


Ko4< 


AV 


A. 


(S)    q  =   KoGo    (1+  Git  +    «s<  )£  - (KoGoH F)  qQ 


where   KoGo  =  Ko    (pt  Zao) 


2- 


P  +    2  (/»,„)  cj*  (WgJ0 

select:    G^      ,    G^       ana  Hp  such   that   substitution 
into      (2)    yields    (1) 

let:    G^  =    A  Za  ;      Gf  =4&   (P  J:  Zao^  Za   \ ; 

(p  f Zao)  Ko  ( p  f Zao) 

„        2AtfHin)p-   A(W§) 
Ko    (p   +  Zao) 

then:      (p2-h  2  ("M|p  +  W§)0      q0=Ko(p+Zao) 
fl+^Za    +AK    (yo  f-  Zao  f  AZa)  )   <f 
(      (°*  Zd^    Ko    (p  +  Zao)  j 

-    {2  A(fsm)    p  4  ,A(W2)}    qQ 
or :  (p2  4-    2[(  f  Wn )  o  f  4  ( /ton )]     p  *  ( W2 ) 0+ A( W§ ) 


(Ko+AK)  (p  f  Zaot^Za)  £ 


which  is  identically  equal  to  equation  (1) 
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Hence,  to  compensate  the  invariant  system,  the 
system  can  be  treated  as  though  the  inner  loops  are  opened. 

Compensation  of  the  Invariant  System 

To  compensate  the  airframe  in  its  invariant  state, 
a  simplified  system  was  employed.   The  servomotor  and  ele- 
vator actuator  dynamics  are  reduced  and  the  dynamics  of 
the  rate  gyro  in  the  feedback  loop  were  ignored.   There 
is  no  loss  of  effectiveness  caused  by  making  these  assump- 
tions.  Since  one  of  the  objects  is  to  compensate  entirely 
for  the  servo  and  actuator  dynamics,  and  the  dynamics  of 
rate  gyros  currently  employed  in  practice  are  such  that 
they  are  of  small  consequence  when  compared  to  the  dynamics 
of  the  other  elements  in  the  system. 


A.   Simplified  System: 


Servo  &   Airframe 
Actuator 


yU-GJ+f 


<?- 


/ 

1? 


Q 
70 


f*6 


Open  Loop  Performance  Equation 

It  is  desired  to  increase  the  damping  ratio  to  a 
value  of  approximately  .7  since  this  damping  ratio  produces 

desirable  dynamic  performance.  It  is  further  desired  to 

2  9 

increase  Wn  to  a  value  much  greater  than  (W^Jo,  since  such 

an  increase  will  nullify,  to  a  large  extent,  the  effect 
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of  parameter  variations  which  occur  when  the  airframe 
is  not  in  its  normal  operating  state,  furthermore,  a  high 
effective  natural  frequency  of  the  seconu  order  system 
ensures  a  "flat"  response  of  this  component  over  a  wide 
frequency  band. 


2« 


f 


{f,^z(fajj.P  *KU  %  =  K  te&J  {&  h  -  ^»  ft  j 


let: 


^^  -     t     (p  *-  2*?  )   t±L     ;    or   the   equivalen 


t   thereof. 


then:  { f  +  [z(faJK)0+   X./tr*]/*    +  (Wf)6  +  <>K*llU  * 

where:  ^-  («/„>;e  *   *0  #,  ^       j    Z  f  »>*  =    z(f^)0  h  X*  Kt 

Choose   values    of  ft  f  ib  such   that  U%  >>^|and  T  ==■  .7 


B.  Determination  of  KcGc: 


let:  ft  C     —  /&fi    I E  t     c  I  ;      or  the   equivalent    thereof. 

Then    the    system   is    effectively   reduced    to: 


'Wl  4- 


* 


if 


where : 


J 


23 


A  frequency  response  analysis  was  made  on  the  re- 
sultant system  employing  a  damping  ratio  of  .7,  an  effect- 
ive natural  frequency  of  20  rad/sec,  a  loop  gain  of  25 
and  a  value  of  Zc  =•  8  rad/sec.   The  open  and  closed  loop 
frequency  responses  are  shown  in  the  chart  below: 

From  the  frequency  plots: 


w 

Open  Loop 

Closed 

Loop 

A#          $ 

M 

N 

.2 

1             2 

.965 

0 

.5 

1             3 

.965 

0 

1.26 

1             6 

.965 

0 

3.10 

1           10 

.965 

0 

8.0 

1.35    12 

.966 

0 

13.0 

1.55      4 

.966 

0 

20.0 

1.6    -20 

.976 

0 

50.0 

.8    -64 

.976 

-2.5 

where:  W   rad/sec.  ,<jf  =  phase  angle ■  J>&-  ;    M  =\~   ^N—  phase  angle  ■$?- 


AR_ 


«% 


)y\. 


?*• 


qo 

Kqm 


Effect  of  aircraft/flight  path  parameter  variations 
on  compensated  system 


h. 


nval 

ent   I 

unct: 

Lona 
» 

l  diai 

jra 

^o 

)        > 

) 

$H 
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As  indicated  on  the  first  page  of  this  Chapter, 
Gr ,  ,G£Z    tiz  Hp~  are  fictitious  functions  usea  to  represent 
the  effects  of  parameter  variations.  The  existence  of  such 
variations  is  equivalent  to  closing  the  loops  containing 
these  functions.   to/ith  these  loops  "closed"  the  following 
analysis  of  the  system  is  made. 

Performance  Function: 

Unity  feedback  loop  open; 

The  frequency  response  of  the  first  bracketed 
term  is  the  open  loop  frequency  response  of  the  compen- 
sated invariant  system.  The  frequency  response  of  the 
second  term  is  determined  by  substituting  the  functions 
for  Gr£j  &  Cjgt  previously  derived,  and  the  resultant  func- 
tion becomes:  K(  (^n/P  *~1-) 

(*k,P+t)  .  ..  . 

where:    /-    =     (Ko  fAKJiZdo  +&2a,) 

The    frequency   response    of   the    third    term   is 
similarly   determined   ana  the    resultant    function   of 

i 


Q+K  $ty)  becomes : 


where : 

2.  _   2 
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Thus,  to  produce  the  open  loop  frequency  response 
of  the  controlled  system,  it  is  necessary  to  determine: 
qo  (  jw)  =  ART/<&- 

where  :   ART  =  log-1  {  log  AR(KoKcGi)G,  f  log  AR(  i  4.  Gl, ,  +  G$-  )  - 
log  AR(l+-KoGJHF),) 

<&r=    ^(KoKcGDGM  +  ^(lfG?1fG^  )  -  M  1+ KoG'HF} 

Using  this  method  of  variations  representation; 
frequency  analysis  plotting  and  transient  response,  by 
analog  simulation,  were  made  for  the  compensated  airframe 
in  seven  different  flight  conditions. 

The  flight  conditions  tested  were  chosen  as  a  re- 
sult of  the  desire  to  impose  the  most  severe  tests  of  the 
effectiveness  of  this  technique  of  compensation.  Limits 
of  altitude,  airspeed,  and  loading  were  determined  from 
current  literature,  and  from  operational  Naval  aviators, 
in  an  effort  to  maintain  the  conditions  within  practical 
limits. 

Definition  of  Flight  Conditions 

Condition  A  represents  the  aircraft  operating  in 
its  invariant  (normal)  state.  Condition  B  represents  the 
aircraft  operating  at  high  altitude,  at  high  speed,  and 
with  payload  expended.   Condition  C  shows  the  aircraft  mak- 
ing an  approach  for  landing,  with  payload  expended.  Condition 
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D  simulates  the  aircraft  at  a  "medium"  altitude  and  at 
reduced  speed,  while  fully  loaded.   Condition  E  represents 
the  aircraft  at  high  altitude,  reduced  speed,  and  payload 
expended.   The  F  condition  shows  the  vehicle  at  the  same 
altitude  and  speea  conaitions  as  E,  but  fully  loaded.  Fin- 
ally, condition  G  simulates  the  aircraft  in  takeoff  condi- 
tion, fully  loaded. 

The  resultant  transient  and  frequency  responses  of 
the  aircraft  in  the  seven  test  states,  both  in  a  compensated 
and  uncompensated  control  state,  for  comparison  purposes, 
are  included  as  Appendices  A  and  C. 
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CHAPTER  5 

MODEL  DESIGN,  VARIABLE  GAIN  AMPLIFIER 
AND  OVERALL  SYSfEM  PERFORMANCE 


Model  Design 

As  indicated  in  Chapter  2,  the  system  design  ob- 
jective is  to  cause  the  system  output  response  to  approx- 
imate as  closely  as  possible  the  response  of  the  model 
to  various  inputs.   The  model  is  the  physical  embodiment 
of  the  design  specifications  of  a  system,  and  a  change 
In  the  model  is  exactly  equivalent  to  a  change  In  system 
specifications.   The  first  consideration,  then,  is  to 
decide  upon  the  system  specifications,  which  must  be 
compatible  with  performance  capabilities  of  the  aircraft. 

The  model  can  then  be  designed  independently  of  the  de- 

(3) 
sign  of  the  control  system. 

Since  the  theoretical  air  frame   was  chosen  with 
no  specific  mission  in  mind,  the  model  reference  system 
for  this  study  was  arbitrarily  selected  to  be  a  second 
order  system  characterized  by  an  undamped  natural  fre- 
quency of  3.5  rad/sec  and  a  damping  ratio  of  .7. 
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Variable  Gain  Amplifier 

It  is  shown  in  Appendix  A  that  the  compensated  con- 
trolled system  response  has  little  phase  shift  from  the 
controlled  system  input,  for  the  different  conditions 
tested  over  a  wiae  range  of  frequencies.   However,  the 
system  does  attenuate  the  input  to  varying  degrees  for 
the  seven  test  condition. 

In  order  that  the  requirement,  that  the  system 
output  response  duplicates  as  closely  as  possible  the 
model  response,  be  fulfilled,  it  is  required  that  an  amp- 
lifier be  inserted  between  the  model  and  the  controlled 
system.   However,  since  the  system  attenuates  its  input 
by  varying  amounts,  depending  upon  the  flight  condition, 
it  is  necessary  that  the  amplifier  be  a  variable  gain  de- 
vice sensitive  to  the  flight  condition.   The  performance 
criterion  of  this  device  is  such  that: 

qp  »  , 

3m 

Employing  the  following  development;   an  acceptable 

performance,  within  the  above  criterion,  of  the  variable 

gain  amplifier  is  obtained: 

Define:  g0  =  ^  ;   gx=  ^  J   Kvar=  gain  of  variable 
yMo  /**  amplifier 

Where  o  subscript  refers  to  the  invariant  state  of 

the  parameter.   The  x  subscript  refers  to  the  general  state 
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of  the  parameter,  and  V,0   and/4  are  as  derined  in  Chapter 


3. 


Required:    Kvar 

or:    Kvar  _ 


25      gx/g 


X/  &o 


1+25    gx/go 
1+25    gx/go 


=  1 


25  gx/ 


X/  SO 


As  indicated  above,  the  use  of  a  variable  gain 
amplifier  having  this  performance  enables  the  controlled 
system  ou,tput  response  to  duplicate  the  model  output  res- 
ponse. 

Overall  System  Performance 

V»ith  the  reference  system  moael  and  variable  gain 
device  established,  a  frequency  response  plot  and  an  analog 
computer  transient  response  were  made  for  the  overall  sys- 
tem.  The  results  of  these  tests  are  included  herein  as 
Appendix  D. 

It  was  found  that  for  all  frequencies  of  interest 
the  overall  system,  i.e.,  the  model,  variable  gain  amplifier 
and  closed  loop  pitch  rate  system;  can  be  represented  by 
the  following  functional  diagram: 
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Where  Katt  represents  the  controlled  system  atten- 
uation of  the  input  to  the  controlled  system  and  Kyar 

Katt  -  1 

This  representation  was  found  to  be  valid  regardless 

of  the  condition  in  which  the  aircraft  was  operating. 
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CHAPTER   6 
CONCLUSIONS  AND   RECOMMJDATIONS 


Conclusions 

Using  only  the  linear  technique  as  described  here- 
tofore, the  pitch  rate  control  system  can  be  compensated 
to  the  extent  that  the  aircraft  response  follows  the  ref- 
erence system  model  output  with  a  phase  shift  of  no  more 
than  ± 10   for  all  practicable  input  rates.   This  method 
of  compensation  is  not  completely  self-adaptive,  since 
the  compensated  system  attenuates  its  input  (the  refer- 
ence system  model  output)  amplitude  in  varying  degrees 
up  to  25$.   Unless  this  attenuation  is  acceptable,  a  var- 
iable gain  amplifier  should  be  employed  to  eliminate  this 
feature  of  the  compensated  system.  Since  the  variable  gain 
amplifier  is  programmed  so  that  it  is  sensitive  to  the 
flight  condition  of  the  aircraft,  it  is  an  overall  non- 
linear device  from  state  to  state.   But,  with  the  air- 
craft in  a  given  state,  the  amplifier  performs  as  a 
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linear  element.  Utilizing  the  combination  of  linear 
compensation  and  a  programmed  variable  gain  amplifier, 
a  pitch  rate  control  system  having  essentially  a  "flat" 
frequency  response,  i.e.,  ±10   phase  shift  and  an  ampli- 
tude ratio  of  1  =t  1%,    is  feasible. 

As  indicated  in  the  scope,  the  basic  objective 
of  this  method  of  compensation  is  the  attainment  of  a 
flat  response.  In  other  words,  the  system  was  designed 
with  frequency  response  characteristics  as  the  principal 
consideration,  i'his  system  when  so  considered  became 
self  adaptive,  to  a  great  degree  of  accuracy. 

Of  secondary  consideration  in  this  design,  but 
of  great  significance  in  actual  practice,  are  the  transient 
characteristics  of  the  compensated  system.  Ideally,  the 
transients  of  the  compensated  system  would  be  of  small 
significance  relative  to  the  reference  system  model  tran- 
sients, or  at  worst,  no  greater  than  the  model  transients. 
Under  this  ideal  arrangement,  the  response  time  of  the 
overall  system  used  in  this  paper  is  the  response  time 
of  the  second  order  model,  i.e.,   .85  seconds  to  come 
within  5%   of  correspondence.  The  analog  computer  tests 
of  the  overall  system,  Appendix  D,  shows  that  this  sys- 
tem does  approach  the  ideal  transient  performance.  How- 
ever, the  computer  results  are  misleading  as  a  result  of 
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mechanization  of  the  problem,  in  that  the  computer 
analysis  was  made  Tor  the  general  case  and  conditions 0 
The  general  case  was  then  extrapolated  by  changing  of 
dial  settings  only,  which  did  not  allow  the  required 
flexibility  for  proper  solution  to  flight  condition  F„ 

The  step  response  equations  representing  the 
air  frame  dynamics  in  its  various  operating  states  as 
tabulated  in  Appendix  C,  show  that  with  the  aircraft 
in  flight  conditions  A,  B,  and  C,  the  transients  are 
of  no  consequence  relative  to  those  of  the  modelo  Con- 
ditions D,  E,  and  G  yield  transients,  which  are  of  the 
same  order  of  magnitude  as  those  of  the  modelo  With  the 
aircraft  in  condition  F  the  response  time  of  the  control 
system  is  noticeably  greater  than  that  of  the  reference 
modelo  Aaaitional  study  of  performance  equations  revealed 
that  the  minimum  value  of  the  quantity  -%j£  required  to 
produce  transients  which  are  iaeal,  within  the  previous 
definition,  is  2030    All  of  the  flight  conditions  tested, 
except  P,  have  values  of  this  quantity  which  are  greater 
than  this  minimum0 

Exclusive  of  the  above  limitation,  the  method  of 
invariant  (normal  state)  parameter  compensation  is  effect- 
ive in  obtaining  self  adaptation  in  the  pitch  rate  control 
system.  It  is  reasonable  to  presume  that  this  same  technique 
will  be  valid  and  useful  in  obtaining  the  desirable  level 
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of  self-adaptability,  in  any  control  system,, 

Recommenaations 

It  is  recommended  that; 

lo    The  system  as  developed  in  this  study  be  mechan- 
ized and  tested  further,.   In  the  mechanization  of  this 
system,  it  is  envisioned  that  a  rate  gyrop  such  as  the 
MIT  10^  rate  gyro,  with  a  static  sensitivity  of  1  mv/mr/sec 
be  employed,.  Such  a  device  ensures  the  attainment  of  gain 
requirements  without  saturation  problems,.  To  attain  certain 
of  the  performance  functions  employed  in  the  development 
in  Chapter  4?  it  may  be  necessary  to  utilize  high  gain 
devices  fed  back  upon  themselves  through  devices  having 
transfer  functions  which  are  the  reciprocals  of  those 
desired,, 

2o    As  an  intermediate  testing  and  analysis  procedure, 
it  is  evident  that  another  type  of  analog  simulation  of 
the  linear  adaptive  system  should  be  made,  using  basic 
mechanisms  and  associated  electronics,,  The  characteristics 
of  these  devices  should  closely  approximate  those  specified 
in  recommendation  #1„ 

further,  some  of  the  amplification  factors  and 
controls  may  be  simulated  by  simple  hydraulic  valves  and 
connectionso 

This  "breadboard"  type  of  study  would  satisfy  a 
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basic  design  goal;  that  of  using  available  components  that 
are  already  capable  of  design  and/or  manufacture „  These 
devices,  once  proven,  can  be  extrapolated  and  modified 
as  designs  for  the  specific   devices  to  be  used  in  the 
control  system  desiredo 

3<,    Additional  study  be  made  on  this  technique  of 
compensation  with  the  objective  of  obtaining  a  compromise 
system  which  yields  an  acceptable  frequency  response  and 
an  acceptable  transient  performance,  i0eos  near  ideals 
for  all  operating  states  of  the  aircrafto 
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APPENDIX  A 

Frequency  Response  of  the  Pitch  Rate  Systemo 

1»  Definition  of  the  Basic  Plight  Conditions 0 

20  Polar  Diagrams  of  the  Compensated  System0 

3o  Polar  Diagrams  of  the  Uncompensated  System; 
Unity  Feedback  onlyQ 
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1.   The  Basic  Plight  Conditions  Defined 


Condi 

tion 

Altitude  (ft) 
55,000 
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S3  (  Slll£ 

J5l 

A 
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a 
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The    GAP/R  Analog  Computer?    General 
Information  &  Techniques » 
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General  Information 

The  GoAo  Philbrick  Researches  Incorporated,  Elec- 
tronic, high  speed  analog  computer  was  used  to  verify 
the  transient  response  characteristics  of  the  uncompensated 
aircraft  pitch  rate  system  with  unity  feedback,  the  com- 
pensated system,  and  the  idealized  model  transient  response0 

It  was  decided  to  test  the  pitch  rate  system,  un- 
compensated and  with  unity  feedback  only,  since  this  is 
an  obvious  first  approximation  to  self  adaptability,  and 
then  compare  this  with  the  response  of  the  final  compen- 
sated unity  feedback  system.. 

The  Computer  used  is  one  of  the  earlier  models  of 
the  GAP/R  unito   The  idea  involved  is  that  of  a  combination 
of  "black  boxes"  or  modules,  that  perform  specific  functions 
within  themselves,  such  as  multiplication,  addition,   inte- 
gration, differentiation,  squaring,  ana  simulation  of  non- 
linear functionso   Further,  these  capabilities  may  be  ex- 
tended to  include  complete  first  and  second  order  servo 
systems c  The  modules  were  designed  to  operate  In  "fast" 
time,  so  that  results  may  be  seen  on  any  standard  oscil- 
loscope trace,  from  which  permanent  photographic  records 
can  be  made „  In  the  modules,  herein  after  referred  to  as 
boxes,  that  are  concerned  with  time  functions,  there  are 
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characteristic  time  constants  (TQ)  settings  of  4  and 
04  raillisecondSo   Modifications  were  made  at  MIT  that 
allowed  "real"  time  operations,  by  changing  the  value 
of  the  condensers  in  the  box  that  is  concerned,,  so  that 
4ms,  04ms  and  1  second,   are  available  for  T0  settingSo 

The  boxes  actually  used  in  the  computer  simula- 
tion are  described  more  fully  in  the  latter  section  of 
this  appendix<> 

The  nominal  voltage  range  for  all  variables  is 
100  volts,  that  is,  minus  50  to  plus  50o   The  manufact- 
urer specifies  that  electronic  variations  affect  the 
sensitivities  by  less  than  \%a      Each  box  is  made  up  of 
one  or  more  DC  minimum  drift  amplifiers  that  have  been 
modified  by  chopper  stabilization;,  The  drift  problem  was 
negligible,  forafter  initial  alignments  checks  on  the 
less  complicated  units  showed  very  little  arift  when 
noted  dailyc   However,  the  more  complicated  systems, 
such  as  the  first  order  lead/lag  networks^,  and  the  sec- 
ond order  generating  system,  which  involves  eight  ampli- 
fiers in  the  units,  required  several  checks  in  the  course 
of  solution0 

As  will  be  seen  in  the  attendant  circuit  diagrams, 
the  operational  amplifier  potentiometers,  which  regulate 
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the  magnitude  of  multiplication  factors,  are  of  the 
proportional  type,  i0eos  one  portion  in  the  feedback 
loop  of  the  bQ    amplifier,  and  the  other  is  in  the  feed 
forward  loop0   -^hese  are  normally  set  so  that  clockwise 
rotation  of  the  potentiometer  dial  puts  more  resistance 
in  the  feedback,  while  reducing  that  in  the  forward  loop, 
thus,  causing  multiplication^,  In  the  K4DY  second  order 
however,  the  dial  A2  is  wound  oppositely,  therefore  div- 
ision is  performed  for  clockwise  dial  rotation   This  is 
necessitated  by  the  form  of  the  standard  second  order 
differential  equationo 

Electronic  clamping  of  the  integrating  boxes  back 
to  zero  potential  is  necessary,  if  the  solution  to  the 
problem  is  not  reached  prior  to  the  characteristic  time 
of  the  integration  units  that  are  involved  in  an  open 
loop  solution,.  If  these  units  are  connected  in  c  losed 
loop,  clamping  is  generally  not  required;   however,,  due 
to  d rifting  DC  levels,  it  is  best  to  clamp  all  integrat- 
ing boxes,  as  a  matter  of  course0 

Calibration  of  the  dial  settings  on  the  multipli- 
cation and  DY  boxes  was  difficult,  since  normal  bridge 
resistance  measurements  would  not  give  satisfactory  re= 
sultSo  Rather,  known  voltages  were  fed  into  the  boxes 
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being  calibrated,  and  the  aial  set  to  give  the  desired 
value  of  output  voltage  „   £his  is  not  too  time  consuming, 
since  other  calibration  raethoas  on  components  require  as 
much  care  and  times  for  accurate  results » 

The  feature  of  high  speed  solution  and  simultaneous 
viewing  is  the  greatest  advantage  of  this  type  of  computer0 
It  is  reaaily  useful  in  simulating  any  control  system,  due 
to  the  many  functions  that  can  be  generated  by  a  series  of 
interconnected  black  boxes,  and  a  quick  solutions  or  quick 
evidence  of  an  improper  solution,  becomes  apparent  almost 
immediately0   Since  so  many  solutions,  to  physical  systems, 
are  in  the  form  of  the  standard  second  order  differential 
equation,  then  the  K4DY  box  inherently  is  ideal  for  a  quick 
approximate  solution0 

Method 

In  setting  up  the  pitch  rate  control  system  on 
the  computer,  there  were  some  modifications  required? 

In  that  the  time  constant  of  the  unit  lead  box  is 
set  from  zero  to  one,  and  it  is  desired  that  all  T0  be 
set  at  o 4ms  (for  ease  in  scaling)  then  the  uncompensated 
open  loop  system  equation  of  the  form: 

qo^KfpT+Dq^ 


p  (2nd  order  equation) 
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is  rewritten  as: 


'  P 


(2nd  order  equation) 

where  1  can  be  simulated  by  an  integration  boXo 
P 

The  same  type  of  problem  arises  from  the  compen- 
sated closed  loop  system  form: 

qo   K(pTl+l)  (pTp-fl)qm  where  Ti)l 
(2nd  order  equation)  / 


and   T2  \  1 


rewrite  as: 


<lo-K(P+.l  )  (pT2+l}Tiqra 
Tl 

(2na  order  equation) 


and  KT-j_  becomes  a  new  K8  9    and  p  can  be  approximated  by 
a  derivative  box0   This  derivative  box,  since  a  perfect 
derivative  is  almost  impossible  to  generate,  causes  a 
time  lag  error  of  1,5  seconds  (problem  time)  in  the  paral- 
lel network  in  which  it  is  operating.   This  leads  to  an 
undesirable  time  lag  in  the  computer  simulation^.  Realizing^ 
however,  that  this  particular  shortcoming  exists;  permits 
its  removal  in  the  final  analysis0  i'his  unwanted  lag  is 
not  apparent  in  the  unit  lead  box» 

The  compensated  system,  has  the  open  loop  form: 
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q0=K(pT+l)  (p  +  Z)qm 

(2nd  order  equation) 

Where  1  approaches  one  of  the  roots  of  the  2nd  order 
T 
equation,  and  Z  approaches  the  other  root.   There  is  high 

gain  in  the  forward  loop,  also.   Then,  to  a  close  approxi- 
mation, the  equation  may  be  written  in  the  final  closed 
loop  form  as: 

q0_K»  (pT-L+Dc^ 

;  K?  includes  K  and  the  T?s 

(pTg  -f-1)  that  are  determined^  as  are 

T-,  and  T2,  according  to  the 
problem  involved, 

Thus,  the  equation  can  be  mechanized  in  the  closed 
loop  by  use  of  a  multiplication  box  in  series  with  a  unit 
lag  and  a  unit  leaa  box,  with  no  feedback  connection  re- 
quired* 

Scaling 

Magnitude  scaling  was  carried  out,  assuming  that 
the  peak  signal  voltage  available  from  the  square  and  sine 
wave  generator  was  13<>5  volts  and  noting  that  50v  was  the 
limiting  factor  of  the  DC  amplifiers,  and  that  the  maximum 
input  signal  to  the  pitch  rate  system  was  from  „1  to  <>3 
radians/secondo   For  scaling  purposes >  the  quantities  E 
and  X  are  introduced  to  the  problem  as  intermediate  quanti- 
ties o 
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E,  the  error,  is  the  difference  between  qm  and  q  <> 
X  is  the  input  to  the  unity  gain  second  order  systenio 

For  simplified  presentation,  the  output  gain  of 
the  system  was  changed  to  show  the  final  results,  in  the 
photographs,  in  the  proper  relative  aspect.  That  is,  the 
output  is  scaled  relative  to  the  input,  i'hus,  if  the  in- 
put is  10  squares,  and  the  output  is  8  squares  in  steady 
stace,  then  the  final  value  is  80^  of  the  input,  and  no 
extra  magnitude  scaling  is  required  of  the  reader0 

i'ime  scaling  for  this  problem  was  set  up  to  be  a 
special  case,  in  that  all  units  involving  time  were  set 
for  T0  equal  to  04ms0  fhus,  1  second  of  "real  time"  on 
the  computer;  in  computer  (problem)  time  is  2500  secondso 

A  dual  beam  oscilloscope  was  used  for  final  pre- 
sentation, in  order  that  both  input  and  output  may  be 
viewed  simultaneously. 

Input  Signals 

response  to  a  step  function,  to  simulate  violent 
maneuvers;  ana  a  ramp  function,  to  simulate  more  gradual 
maneuvers,  was  used.   The  frequency  of  the  square  wave 
step  function  was  30  cps.  In  problem  time,  this  is  42 
seconas/step.  i'he  slope  of  the  negative  input  ramp  is 
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o 12  units  of  input/computer  second.  Each  small  square 
on  the  grid  of  the  photographs  represents  lo67  seconds 
of  problem  time0 

Frequency  Response 

The  scope  of  this  paper  does  not  include  experi- 
mental frequency  analysis.,   However,,  in  preliminary  study 
of  the  problem  it  was  found  that  the  more  complicated 
boxes  have  their  own  frequency  characterist ics0   It  was 
also  noted  that  to  make  proper  response  for  lead,  or  dif- 
ferentiating systems,  the  following  must  be  done«   Connect 
a  system  to  s  imulate  a  lead  lag  transfer  function  made 
up  of  integrators,  potentiometers  and  adders,,  Make  the 
denominator  term  of  the  network  far  enough  out  on  the 
real  axis,  so  as  to  essentially  leave  a  pure  lead  term 
in  the  final  result,  that  gives  the  best  frequency  res- 
ponse„   The  unit  lag  boxes  exhibited  better  frequency  char- 
acteristicso   The  second  order  box  was  accurate  for  the 
frequencies  tested,  and  should  adapt  well  to  measurement 
of  simulated  sinusoidal  gust  type  responses^ 

Circuit  Representation 

A  more  detailed  stuay  of  the  unit  boxes,   their 
functions  and  capabilities,  along  with  a  shorthand  method 
of  representation,  is  presented  in  later  sections  for  the 
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following  reasons: 

lo    Since  most  present-aay  manufacturers  of  analog 
computers  have  all  but  abandoned  the  modular  technique; 
the  internal  circuitry  is  shown  for  comparison  with  present 
simplified  "plug  in"  R  and  C  component  methods0 

2o    To  demonstrate,  in  block  diagram  form,  the  com- 
puter final  circuit  connections;   to  simulate  the  pitch 
rate  system,  model,  and  associated  compensation  deviceso 
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APPENDIX  B-2 


Circuit  and  Block  Symbols  of  the 
GAP/R  Components o 
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KZ-A  Adding  Component 


Equation:   y  *  x±+  xg  +  . .  .X 


Block  Symbol: 


4 


>.+y 


Note: for  simplicity, the 
JIC  and  NEMA  symbol  for  control 
resistance  is  utilized: 


DESCRIPTION:   The  "Adder"  computes  the  instantaneous  sum 
of  four  or  fewer  input  signals.  By  interconnecting  Adders 
in  tandem,  any  number  of  signals  may  be  combined  in  sums 
and  differences.  V»ith  N  such  Components,  3N  -t-1  signals 
may  be  added.   Each  Adder  features  a  calibrated  dial  where- 
by a  steady  voltage  may  be  aaded  in,  amounting  at  the 
extremes  to  plus  or  minus  10  volts,  or  20$  of  the  maximum 
signal  excursion.   This  is  convenient  for  manual  addition 
of  a  constant  term  in  an  equation. 

Signals  receive  neither  gain  nor  attenuation  through 

the  Adding  Component. 

57 


K3-C  Coefficient  Component 


Equation: 


y=  ±Cx 


100  -,-0 


-y 


Block  Symbol: 


x >  c     y 

— ^ 


DESCRIPTION:    An  input  signal  to  this  Component  is 
automatically  multiplied  by  a  numerical  factor,  then 
carried  directly  and  inversely  to  the  two  outputs.  The 
factor  may  be  set  on  a  specially  calibrated  nonlinear 
dial,  extending  from  0  to  100.   Technically,  gains  from 
a  precise  zero  to  infinity  are  within  range,  but,  the 
best  linear  range  is  0.1  to  10. 
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K3-J  Integrating  Component 


Block  Symbol: 


DESCRIPTION:   This  Component,  the  "Integrator",  computes 
an  integral  with  respect  to  time  of  the  voltage  supplied 
to  the  lower  input.  The  time-factor  T   is  0.4  millisecond. 
For  normal  computing,  in  closed  loops,  the  right-hand 
switch  position, is  chosen  for  maximum  sensitivity;  the 
left-hand  switch  position  introduces  a  stabilizing  net- 
work and  provides  for  open-cycle  operating  when  desired, 
as  for  the  purpose  of  display  alone. 

The  integrated  outputs  are  artificially  and  auto- 
matically returned  to  zero  by  means  for  a  "clamping"  sig- 
nal supplied  to  the  upper  input. 
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K3-D  Differentiating  Component 


Eauation: 


ty=  T0dx  =  T0px 
dt 


+  y 


Block  Symbol: 


DESCRIPTION:    Operationally,  the  "differentiator" 
performs  inversely  to  the  Integrator.  It  computes 
time  derivative  of  the  input  signal.  The  time-factor 
T0  is  0.4  milli-second.  With  the  switch  at  the  left- 
hand  setting,  as  in  this  case,  a  lag  is  introduced 
having  a  time  constant  of  about  50  microseconds; 
about  l/lO  the  fixed  time-factor  of  the  Component. 
This  setting  is  an  approximate  derivative,  which 
did  affect  accurate  representation  showing  as  spurious 
lag. 
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Ko-E  Augmenting  Differentiator 


Equation:  +  y»  x  +  Tdx  s  (l+Tp)x 

dt 


Block  Symbol: 


DESCRIPTION:   This  Component  is  related  in  a  simple 
way  to  the  K3-D  Differentiator,  since  it  passes  the 
input  signal  and  adds  to  it  the  first  time  derivative 
of  that  signal.   The  time-factor  T,  the  derivative 
sensitivity,  is  continuously  adjustable  from  zero  up 
to  0.4  milli-second  on  a  linear  percentage  scale. 
The  operator  embodied  in  inverse  to  that  of  the  K3-L 
Unit-lag  Component. 

Pbr  certain  purposes  this  Component  may  replace 
a  combination  of  K3-A,  K3-D,  and  K3-C.   It  represents 
the  introduction  of  a  derivative  response. 
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K5-L  Unit-Lag  Component 


Equation:    y  ♦  Tdy_  =(l  +  Tp)y=±x 

dt 


Block  Symbol: 


>  ♦ 


DESCRIPTION:   The  "Unit-lag"  performs  a  first-order 

lagging  operation.  By  means  of  its  0-100  dial,  the 

time  constant  T  may  be  adjusted  from  zero  up  to  a 

maximum  of  0.4  millisecond.   The  final,  or  long-term, 

sensitivity  is  unity.   This  type  of  operation  or 

dynamic  characteristic  may  be  obtained  by  a  simple 

loop  involving  a  K3-A,  a  K3-J,  and  a  K3-C,  which 

loop  also  affords  several  other  characteristics. 

However,  the  K3-L  Component  provides  a  compact  and 

convenient  version  of  this  frequently  recurring 

operator. 
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K4DY  Dynamic  Component 


Equation:    A2T§  d^y  +  AXT0  dy_  +  A0y  =  x(t) 


Block  Symbol: 


DESCRIPTIOM:   This  box  combines  the  functions  of  eight 
assorted  A,C,  and  J  boxes  into  one  simplified  unit.   The 
standard  second  order  equation  satisfied  by  this  unit, 
is  an  intrinsic  part  of  most  actual  physical  systems. 
TQ  is  equal  to  .4  milliseconds,  and  Ai,Ag,A0  are  directly 
calibrated  on  0-100  dials  of  the  Coefficient  type.  Dial 
A2  is  wound  inversely,  as  previously  described. 
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APPENDIX  B-3 

GAP/R  Block  Connection  Diagrams  and  Dial 
Settings,  for  the  Pitch  Rate  System  Simulation,, 
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FLIflHT  CONDITION  A 
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Computer  Block  Diagrams  for  the  Compensated  Pitch  Rate 
System, During  the  Various  Flight  Conditions. 
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Computer  Block  Diagrams  for  Uncompensated  Unity  Feedback 
during  the  Various  Flight  Conditions, Pitch  Rate  System. 
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Computer  Diagram  and  Dial  Settings  for   the  Model 
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Unit  Dial  Settings  for  Compensated  System, 


\Uni 
Flight^^ 
Condition^ 

t 

VC1 

G2 

G3 

E 

L    A2 

Al 

Ao 

A 

,185 

,500 

1 

ol25 

ol 

2,8 

4 

B 

,936 

ol26 

,123 

C 

,999 

ol25 

,125 

D 

8ol4 

4,37 

ol24 

5 

4,56 

lo03 

E 

8,43 

4o82 

,125 

5 

4,43 

1,03 

P 

80  02 

2,32 

,125 

5 

4,71 

1,44 

G 

7,98 

6,15 

,125 

5 

4,73 

lo43 

a 


K4DY  and  K3D  set  for  To  =  ,0004  sec. 


b0   A2,  A-,  and  A0  are  the  dial  settings  of  the 
K4DY  2nd  order  equation  module. 
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Unit  Dial  Settings  for  the  Uncompensated 
System,,   (with  unity  feedback) 


sUnit 

Flighl 

Condition 


E 


A 

B 

C 

D 

E 

F 

G 


.925  .925  o69 

loll  1.00  o073  o094 

,925  0925  ol67 

loll  loOO  o018  o034 

1011  loOO  o031  .054 

loll  loOO  oOl  e025 

1.11  loOO  o02  .025 


935  2.5   5.25  lo 21 
.500  o300   0092 

325   2.5  15ol  2,32 
.500  0214  o033 
.500  o230  .05 
.500  .115  .025 
.500  .30   .026 


a.   K4DY  and  K3J  are  set  for  T0  =  .0004  sec. 

bo   Ag,  A-j_,  and  AQ  are  the  dial  settings  of  the 
K4DY   2nd  order  equation  module. 
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APPENDIX  C 

Transient  Reponse   of  the  Pitch  Rate  System., 
lo   General  information  and  definitions0 
2,   Step  input,  simulating  violent  maneuver, 
30   Ramp  input,  simulating  gradual  maneuver,. 
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GENERAL  INFORMATION  AND  DEFINITIONS 

lo   Theoretical  step  response  of  the  compensated  pitch 
rate  system,  determined  by  inverse  transform  methodso 

Flight  Condition  Step  Response 

A  qQ^  o965  +o0222  e  "8o2t 

B  q0  =  0945  +  o0245  e  ~8°155t 

C  qo=o990-  o00497  e  "8o04t 

D  q0_.768-|-  o092  e   °  toU^+o0545e 

E  qQ=  084  +.0357  e  8o3t  -K073  e  =  o62t 

F  q0^  06-|-0268e~°425t-K0552e=8o5t 

G  q0=.  o  761-j-  O0587e  ~o83t-+ O067e°8°  64t 

20   The  uncompensated,  unity  feedback,  system  reduces  to 
a  closed  loop  type  1  servo  systenu   The  final  value  theorem 
predicates  that  the  final  amplitude  ratio  of  output  to  input 
is  equal  to  1„ 

3o   Each  small  block  on  the  photographs  is  equal  to  lc67 
seconds  of  problem  time0 

4o   For  the  ramp  input  function;  the  output  ramp  is  ex- 
trapolatea,  ana  the  "dynamic  error"  is  estimatedo 
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Dynamic  Error  (DE)  is  defined  as  the  difference 
between  the  output  response  level  at  fe,  and  the  forcing 
function  point  for  the  same  instant: 
or  in  equation  form: 

(DE>*i=<l(out}t>l  -  *(out)  qs  tx 
where  qs =  quasi-static,  the  estimate  of  the  output 

level  rate,  for  quasi-static  conditions  • 

(DE)  is  then  nondimensionalized  as    /INPUT 

5,    R>r  the  step  input  function,  the  "response  time"  is 
determined,, 

Response  time  (RT)  is  defined  to  be  the  time  required 
for  the  output  response  to  change  from  an  initial  arbitrary- 
initial  dynamic  error  to  within  five  percent  of  the  forced 
output  response,  or  to  within  9b%   of  its  steady  state  value0 
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(A)  Compensated 


RT   .835  seconds 


DE/Input   .05 


(A)  Uncompensated 


RT   5.85  seconds 


DE/Input   .1 
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(B)    Compensated 


RT      .167   seconds 


DE/Input      .1 


(B)   Uncompensated 


RT     10.85   seconds 


DE/Input      .561 
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(C)  Compensated 


RT   .354  seconds 


DE/Input   .1 


(C)  Uncompensated 


RT  10.85  seconds 


DE/Input   .561 
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(D)  Compensated 


RT  5.0  seconds 


DE/Input   .474 


(D)  Uncompensated 


RT  15.4  seconds 


DE/Input  .5 
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(E)    Compensated 


RT     7.51   seconds 


DE/Input      .5 


(E)   Uncompensated 


RT     10.85   seconds 


DE/Input      .208 
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(F)    Compensated 


RT     4.58   seconds 


DE/Input      .465 


(F)   Uncompensated 


RT     15.08   seconds 


DE/Input    .615 
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(G)  Compensated 


RT   6.51  seconds 


DE/Input   .25 


(G)  Uncompensated 


RT  17.5  seconds 


DE/Input   .80 
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APPENDIX  D 

Transient  and  Frequency  Response  of 
the  Combined  Model  and  Pitch  Rate 
Systemo 
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The    performance    function   of  the   model   in 
closed  loop   is: 


%i     wn  1  in 


p2f  2fWnp  +  wg 


where:      Wn=3o5    radians/second 
/  =  .? 
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Combined  System 


RT   1.1  seconds 


DE/lnput    0 


*A'*=\&\ 


1.0 


APPENDIX  E 


BIBLIOGRAPHY 


83 


1„    Rom,  Meeker  and  Lund,  Self -Optimizing  Aircraft 

Pitch  Control;  Report  T-119,  Instrumentation  Lab- 
oratory,  MIT;  MoS*  Thesis;  May,  19570 

2o    Osder,  S.S-  and  Hutchinson,  IoNa,  Pinal  Technical 
Report  Feasibility  Study  Automatic  Optimizing  Sta- 
bilization System,  Part  1;  WADC  Technical  Report 
58-243  Part  1,  ASTIA  Document  AD  155576;  June,  1958Q 

3o    Whitaker,  Yamron  and  Kezar,  Design  of  Model-Reference 
Adaptive  Control  Systems  for  Aircraft;  Report  R-164, 
Instrumentation  Laboratory,  MIT;  September,  19583 

40    Nonweiler,  To,  Descent  from  Satellite  Orbits  Using 
Aerodynamic  Braking;  from  Realities  of  Space  Travel 
editor,  L.J*  Carter;  McGraw-Hill;  New  York,  19570 

50    Schmidt,  SQFa,  and  Triplett,  ViuCo.,  Use  of  Non° 
linearities  to  Compensate  for  the  Effects  of  a 
Rate  Limited  Servo  on  the  Response  of  an  Automat- 
ically Controlled  Aircraft;   Ames  Aeronautical  Lab- 
oratory, Moffat  Field,  Califprnia;   NACA  Technical 
Note  3387;  January,  1955o 

7o    Flugge-Lotz,  K-.  ,  and  Taylor,  C.F.,  Investigation 

of  a  Nonlinear  Control  System;  Stanford  University; 
NACA  Technical  Note  3826,  April,  1956'0 


84 


80    Lowes  A.,  and  Danois,  Mo,  Analysis  and  Synthesis 


of  a  Linear,  Self-Adaptive,  Stability  Augmentation 
System;   Aerophysics  Research  Report  FZA-271,  Convair; 
January,  1959 „ 

9,    Locke,  A.S.,  Guidance,  D«  Van  iMostrand  Company, 
New  York,  1955o 

10 o   McLane  and  Van  Meter,  A  Study  to  Determine  the  Auto- 
matic  flight  Control  Configuration  to  Provide  a 
Stability  Augmentation  Capability  for  a  High  Perform- 


ance Supersonic  Aircraft;   Minneapolis-Honeywell; 
WADC  Technical  Report  57-349,  Parts  1,  2  and  3o 
ASTIA  Documents  AD  156146,  156145,  157753o 

llo   Larrabee,  EaEo ,  Unpublished  notes  for  a  class  in 
Advanced  Airplane  Stability  and  Control;  Depart- 
ment of  Aeronautics  and  Astronautics,  MIT;  September, 
1958c 

12o   Li,  YaTo.,  The  General  Research  Study  of  Dynamic 

Response  Self-Adaptive  Control;  Technical  Proposal; 
Department  of  Aeronautics  and  Astronautics^,  MIT; 
December,  1958o 

13o   Wheeler,  RaCaHa,  Unpublished  notes  for  a  class  in 

Servomechanisms;   Department  of  Electrical  Engineer- 
ing, UaS*  Naval  Postgraduate  School;  October*  19570 

85 


14o   Thaler,  G-.  F. ,  Unpublished  notes  for  a  class  In 

Advanced  Servomechanisms ;  Department  of  Electrical 
Engineering,  U.S.  Naval  Postgraduate  School; 
January,  1958 . 

150   Stein,  Y*»A«,  Unpublished  notes  for  a  class  in 

Electrical  Transients  ana  Analog  Computers;  De- 
partment of  Electrical  Engineering,  U„Sa  Naval 
Postgraduate  School;  August,  19570 

16.   Paynter,  H.M.,  A  Palimpsest  on  the  Electronic 

Analog  Art;  G-aA.  Philbrick  Researches  Inc0,  Boston; 
1955c 


86 


' 


thesA428 

A  study  of  the  linear  self-adaptive  syst 


3  2768  001  89851  3 

DUDLEY  KNOX  LIBRARY 


ll|liiffli^H|ffl 

Html 

l7  trfl I 

II 
IJj 

lit 

'Mjvpf'llllflJ  iiiilntnilHriflKiiiiuliiiiHlriill 

•  ill}     '         i   t[fi}ll     '    '■   }|}i]  '  '■■'  HI]  PlK'tiili 

• 

- 

IliSI™ 

F|||  l|  lif  t{|{  If  ml  IIP  Rlf  iill«i*Ii£i%i 

'<  •  >i  If  Mi  hh  WMnni !  nl'Mi" 

u?]itlt([»ifft»JnilH|fKu|{jri}tf]{](( 


lilt] 


■Hn 


:h( 


[flnjlj 
ijljjjl 

llllll 
j/IIil 

1        i llfSi  If illflBB 

:?  1 1  ||H  f  f  k/I  ff  }«fti*l  II  f  1  If  f  f  fl  flic  r 

ll<    :  \  '  '■■> flwiiltllliH  f':,;i 

1         ISS    f 1 
(lillFllllill  lilliiinlitiiliillilili 

BlIIUHH 

H            P  B  H        1 

IHISSSiiiffi 

I  ml  '  mm 

111 

:    if   I   1/ ml Im 


Hi 


i    » 


n<i!i  illtf! 


11    !f 
■piil:><;ii.< 


